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Ground state properties of the S — 1/2 antiferromagnetic XXZ chain with Fibonacci ex- 
change modulation are studied using the real space renormalization group method for strong 
modulation. The quantum dynamical critical behavior with a new universality class is pre- 
dicted in the isotropic case. Combining our results with the weak coupling renormalization 
group results by Vidal et ai, the ground state phase diagram is obtained. 
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1. Introduction 

The magnetism of quasiperiodic systems has been one 
of the main subject of continual studies since the dis- 
covery of quasicrystals. 1 These quasiperiodic systems 
have an intermediate character between the regular and 
random systems. Among them, the single particle elec- 
tronic states in the quasiperiodic chains which are the 
one dimensional counterpart of the quasicrystals have 
been intensively studied from theoretical viewpoint. 2-5 
Recently, this problem has been attracting renewed in- 
terest following the synthesis of magnetic quasicrystals 
with well-localized magnetic moments. 6,7 The theoreti- 
cal studies of the quantum magnetism in one and two 
dimensional quasiperiodic systems are also started by 
many authors. 8-18 

It is well known that the S — 1/2 XY chain can be 
mapped onto the free spinless fermion chain. The prob- 
lem is, however, not trivial on the quasiperiodic lattice. 
This problem has been extensively studied by Kohmoto 
and coworkers 2-5 by means of the exact renormalization 
group method. It is shown that the ground state of the 
XY chain with Fibonacci exchange modulation is crit- 
ical with finite nonuniversal dynamical exponent. This 
approach was extended to include other types of aperi- 
odicity and anisotropy. 8 It is clarified that the criticality 
of the Fibonacci XY chain emerges from the marginal 
nature of the Fibonacci and other precious mean aperi- 
odicity. For the relevant aperiodicity, more singular be- 
havior with a divergent dynamical exponent is realized 
even for the XY chain. 8 

The investigation of the S — 1/2 Fibonacci Heisen- 
berg chains started quite recently. In the fermionic lan- 
guage, the Ising component of the exchange coupling cor- 
responds to the fermion-fermion interaction, leading to 
the strong correlation effect which is one of the most im- 
portant subject of recent condensed matter physics. In 
the weak modulation regime, Vidal and coworkers 12 ' 13 
have shown that the Fibonacci modulation is relevant 
on the basis of the weak coupling renormalization group 
(WCRG) calculation. The present author carried out the 
DMRG calculation and investigated the scaling proper- 



ties of the low energy spectrum. 9, 10 

In the present work, we investigate the 5 = 1/2 Fi- 
bonacci XXZ chains using the real space renormalization 
group (RSRG) method 19 which is valid for the strong ex- 
change modulation. Combining the RSRG and WCRG 
results, the schematic ground state phase diagram is ob- 
tained. Some of the results for the isotropic case are 
briefly reported in ref. 11. 

This paper is organized as follows. The model Hamilto- 
nian is defined in the next section. The real space renor- 
malization group formulation is given in section 3. In sec- 
tion 4, the solution for the recursion equation is obtained. 
A new universality class is predicted for the isotropic 
case. The ground state phase diagram is presented. The 
last section is devoted to summary and discussion. 

2. Model Hamiltonian 

Our Hamiltonian is given by, 

N-l 



H = 



i=i 



J [ qx QX , qV qV 



A ai S!Si +1 \(l) 



( J a% > 0, on = A or B), 



where Si are the spin 1/2 operators. The exchange cou- 
plings J ai (= Ja or Jb) and anisotropy parameters A ai 
(= Aa or Ab) follow the Fibonacci sequence generated 
by the substitution rule, 



AB, B -> A. 



(2) 



In the following, we consider the initial Hamiltonian 
with Aa = Ab = A. However, the renormalized values 
of Aa and Ab are not equal to each other. Therefore we 
keep the a^-dependence of A Qi in eq. (1). 

3. Real Space Renormalization Group 

If one of the couplings Ja or Jb is much larger than the 
other, we can decimate the spins coupled via the stronger 
exchange coupling and calculate the effective interaction 
between the remaining spins by the perturbation method 
with respect to the weaker coupling. 19 In the following 
subsections, the cases Ja << Jb and Ja >> Jb are 
discussed separately. 
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Fig. 1. The decimation procedure for Ja << Jb- 



Fig. 3. The decimation procedure for Ja » Jb. 



ABAABABAABAABABAABABAABAABABA ■ 
/ \ \ t It I t I t 

^ B7 A' B' A' A' B' A' 



Fig. 2. The real space renormalization scheme of the Fibonacci 
XXZ chain. The letters A and B correspond to the bonds and 
the up and down arrows to the spins which survive decimation. 
For Ja >> Jb the end spin and bond in the parenthesis do not 
appear. 



3.1 Case J A « Jb 

In this case, the Jb bonds are decimated as shown in 
Fig. 1. The spin- 1/2 degrees of freedom survive on the 
sites in the middle of the sequence AA. The two kinds 
of sequences of bonds ABA and ABABA arc allowed be- 
tween two alive spins. There exists one singlet pair in the 
former sequence while two singlet pairs exist in the lat- 
ter sequence. Therefore the effective coupling is weaker 
for the latter case. This decimation process replaces the 
sequence ABABA by A' and ABA sandwiched by two 
As by B' resulting in the sequence B'A'B'A'A'B'A'B'A' 
.. as shown Fig. 2. Except for the B' at the leftmost po- 
sition, this again gives the Fibonacci sequence as proven 
in Appendix. The effective coupling can be calculated by 
the straightforward perturbation theory as, 



Ji = 



A 



r 3 

J A 



An = 



J|(1 + A B ) 2 

Jb(1 + A b ) 

Aj(l + A B ) 2 
4 

A|(1 + A B ) 



(3) 



up to the lowest order in Ja- 



3.2 Case J A » Jb 

In this case, the Ja bonds are decimated as shown in 
Fig. 3. The three spins connected by the successive A 



bonds form a doublet which can be described as a single 
spin with magnitude 1/2. Therefore the same replace- 
ment as the preceding subsection works as well. In this 
case, no B' appears at the leftmost position as shown in 
Fig. 2. Again, two kinds of sequences of bonds are allowed 
between two alive spins, namely B and BAB. There is no 
singlet pair in the former sequence while one singlet pair 
exists in the latter sequence. The effective coupling is 
calculated by the perturbation theory as, 



A = 



Jb 



A' = 



A' = 



4JI 



J A (1 + A A )(8 + A1) 
4J b 



; + a 2 



A 2 



(Aa + \/8TA^) 2 (1 + Aa) 
32 



A B (A A + v /8TA^y 
16 



(4) 



up to the lowest order in Jb- In this case, the ratio of 
two couplings is reversed after the first decimation as, 

Jb 



J4 



<< i. 



(5) 



, B J A (1 + A A ) 

Therefore, after the second decimation, the formulae in 
(3) apply. 
We thus have in general, 



j(n)3 



J<f )2 (l + A B (">) 



2 ' 



J 



(n+1) _ 



J 



(n)2 



fn+1) 



r(«) 



A (n)3 



1 + A 



(n) 



(i+A^y 



(6) 



(7) 



(8) 
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with 



A (n+1) A V B J 



x(l) : J A 

J 2 (1 + A B ) 2 ' 

t(1) j a 

Jb - J B (1 + A B )' 



(9) 



Ai(l + A B ) : 



(1) Ajg + Ag) 

B 2 



for Jb >> Ja and 



J 



(i) 



4J1 



J 



J A (1 + A A )(8 + A A )' 
4J b 



(i) 

B 



+ A 2 ' 



A 



A 2 (a a + 78 + A^) 2 (1 + A a ) 



32 



A B (A A + v/8 + A^)' 
16 



for J A >> Jb where the variables with (") refer to the 
values after n-th decimation. 

4. Solutions of the Recursion Equations 

4-1 XY and isotropic case 

In the XY (A A = A B = 0) and isotropic (A A = Ab = 
1) case, the recursion equations (8) and (9) imply that 
A A and Ab are not renormalized. Therefore the recur- 
sion equations are simplified as, 



t("+1) 

j a 



J 



(")3 



r(n)2' Jb 



(n+1) _ J A 



(n)2 



J 



(n) 



(n > 1) (10) 



with 



J A t(1) _ ^A 
2 ' ~ J B 



J 



J B T (l) _ 
2J ; 



for the XY case and 



r("+l) 
J A 



r(»)3 



for J A << Jb 
for J A >> J B 

r (n)2 



4J 



j£+1) = Ja_ („>1)(H) 



(n)2' J B 



2J 



(n) 



with 



i |, J« = ^ for Ja«Jb 
2J 2 (1) _ 4J B 



9Ja' Jb 9 f ° r Ja>> Jb 

for the isotropic case, 

These recursion equations can be cast into the linear 



recursion equations for hiJ A "^ and lnJ B "^ as follows, 



= MI (n) a for n > 1, (12) 



where 



x (n) 



InJ 
lnj 



(n) 
A 

(n) 



M = 



3 -2 
2 -1 



The constant vector a — (a a, «b) is given by 



a 



in the XY case and 



a = 



1 



ln2. 



in the isotropic case. 

This equation can be solved easily by noting that 



M 



O = 



I + O 



2 -2 
2 -2 



(13) 
(14) 

(15) 
(16) 



and I is the identity matrix. Because O 2 = 0, the follow- 
ing identity holds 

M" = (I + 0) n = I + nO. (17) 

The n-times iteration of eq. (12) gives 

n-2 



k=0 



n-2 



= (I +(n-l)0)X {1) + k °) a 



k=0 



= (I + (n - 1)0) X« 



(n-lU+ ^-^-^ Ola (18) 



Namely, we have 



InJ 



(n) 
A 

r(») 



2n - 1 -2n + 2 \ / InJ 



r(i) 



(n - l) 2 -n 2 + 3n - 2 

n 2 - 3n + 2 -n 2 + An - 3 



(19) 



for n > 1. In the following, we discuss the XY and 
isotropic cases separately. 

4.1.1 XY case 

In the XY case, a = 0. Therefore, the last term of (19) 
vanishes and the solution is given by, 



InJ 
lnj 



(n) 
A 
(n) 



2n - 1 -2n + 2 \ / InJ 



2n - 2 -2n + 3 



lnj 



A 1} 



(20) 
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Namely, 



(«) 



J 



J 



(i) 



(i)' 

A 



2(n-l) 



J 



{n) 



J 



(1) 



(1)\ 2(n-l) 



(1) 



For Jb >> Ja, this reduces to, 



An) _ 



= Ja 



An) 
J B 



Jb 



and for J B << Ja, 



J 



(n) 



J 



(n) _ 



2Ja 

Jb 
2 



Ja 
Jb 



Jb\ 
Ja/ 



2ra 



2n 



2(n-l) 



Jb\ 
Ja) 



2(n-l) 



(21) 
(22) 

(23) 
(24) 

(25) 
(26) 



It follows that the ratio Jp^/J^ 



is invariant under 



rcnormalization after second step of rcnormalization as, 



r("+l) 

[a 

r(«+i) 



J 



(n) 



J 



(n) 



(1) 
A 

(1) 



' AB- 



(27) 



where tab = max ' 7a ) ' Therefore, the aperiodic- 
ity is marginal for the XY case. This is consistent with 
the result of the exact renormalization group method 2-5 
and weak coupling renormalization group method. 12 ' 13 
This also implies that the perturbation approximation 
remains valid if it is in the first step of renormalization. 
As explained in Appendix, our decimation rule 

ABABA -> A, ABA -> B 

is almost equivalent to the usual 3-step deflation process 

ABAAB -» A, ABA -> B 

except for a single B which remains at the left end. It 
should be noted that the length of the 3n-th Fibonacci 
sequence is equal to the Fibonacci number F^, n which 
grows as § Zn for large n where (j> is the golden mean 
(= 1+ 2 )■ This implies that the finite Fibonacci sequence 

of length TV reduces to a single spin pair coupled via 
after n = \nN/ (31n</>) steps of decimations. Therefore, 
the smallest energy scale AE for the finite Fibonacci XY 
chain with length N scales as, 

2 r lnr AB . (28) 



AE 



with z 



31n</> 



This reproduces the result of Kohmoto and coworkers 2 5 
in the limit tab >> 1- 



4-1.2 Isotropic case 

In the isotropic case, a =* (2, l)ln2 and the solution 
is given by, 

2n-l -2n + 2 
2n-2 -2n + 3 

-n 2 + in - 
-n 2 + 4n- 



It should be noted that the last term is most dominant 
in the limit n — > oo. We thus have, 

2(n-l) 





(n) 



J 



J 



(1) 



(1)' 

A 



A 
j(n) 



J 



J 



(1) 



B 



(1)' 
A 



2(n-l) 



For Jb >> Ja, this reduces to, 



7 {n) 



An) 



Ja 



= Jb 



Ja 
Jb 

Ja 
Jb 



2n 



-(n-lY 



^(n+l) 



and for J B << Ja, 



An) 

J A 



An) _ 



8J B 

9 

8J A 
9 



2 



Jb 
Ja 



JbV 
J a) 



It follows that the ratio J^ / J^ 1 



-n(n-l) j (3 ) 

(31) 

(32) 
(33) 

(34) 
(35) 

decreases under 



-n(n+l) 



renormalization as, 

t(«+1) 



1 J 



(«) 



An+1) 



2 An) 
J B 



(n > 1) 



and for the first step, 



(i) 

A 
(1) 



"AB- 



(36) 



(37) 



Therefore the aperiodicity is relevant in the isotropic 
case. This is consistent with the result of the WCRG. 
This also implies that the perturbation approximation 
becomes even more accurate as the renormalization pro- 
ceeds. 

Taking into account that n = ln7V/(31n0) the small- 
est energy scale AE for the finite Fibonacci chain with 
length N scales as, 



AE 



2~ ; 



- exp(-(lnA731n</)) 2 ln2) 
<ln7V with k = In2/(31n0) 2 (38) 



for large enough N, irrespective of the value of tab- This 
type of quantum dynamical critical behavior, which im- 
plies the logarithmic divergence of dynamical exponent, 
is not yet reported in any regular, quasiperiodic and ran- 
dom systems in the field of one-dimensional quantum 
magnetism. 
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This size dependence of the energy scale implies that 
the number of the magnetic excited states ND(AE)dAE 
in the energy range AE ~ AE + dAE can be written 
in the scaling form in terms of the scaled variable x = 

ND(AE)dAE ~ / (x) dx 



Nf iVexp 




exp 



2kAE 



dAE 



(39) 



with a scaling function f(x). Because the density of state 
per site D(AE) should be finite in the thermodynamic 
limit N — > oo, the scaling function f(x) should tend to 
a finite value as x — > oo (N — > oo). Therefore we find, 



D(AE)dAE 



exp 



2kAE 



dAE (40) 



for large enough N. Accordingly, the low temperature 
magnetic specific heat C{T) at temperature T should 
behave as, 

C{T) ~ ~&T N J AED ( AE ) dAE ~ NTD(T) 



N 



2k 



ilnJ- 

K T 



: exp 




(41) 



The magnetic susceptibility x(T) at temperature T 
should be the Curie contribution from the spins alive 
at the energy scale T. The number n s (T) of such spins 
is given by 



n s {T) ~ 27V f D(AE)dAE, 
Jo 



(42) 



because two spins are excited by breaking a single ef- 
fective bond with effective exchange energy less than T. 
Therefore the low temperature magnetic susceptibility \ 
behaves as, 



X(T) 



2N 
4T 



J" 

Jo 



exp 



2nAEj^\n^ E 



dAE 



iVexp(- v /IinT) 



2T 



(43) 



In ref. 9, the present author carried out the DMRG cal- 
culation for the Fibonacci antiferromagnetic Heisenbcrg 
chains. We performed the finite size scaling analysis of 
the lowest energy gap AE based on the assumption that 
it will behave in the same way as the XY chain with rel- 
evant aperiodicity, namely as AE ~ exp(— cN w ). How- 
ever, the present analysis suggests the different behavior 
(38). Therefore we tried to replot the previous data in 
ref. 9 using the scaling (38). Unfortunately, however, the 
fit is very poor. The reason will be understood in the 



following way. 

In the present RSRG analysis, we started the RG 
transformation based on the perturbation theory assum- 
ing that one of the exchange coupling is much larger than 
the other (r\B >> 1)- However, the DMRG calculation 
in ref. 9 is carried out for relatively small tab ■ Therefore 
in the early stage of renormalization, the strong cou- 
pling recursion formulae are not valid. To supplement 
this point, we introduce the phenomenological modified 
recursion formulae 



An+l) 



j(n)3 

-^exp(-a A ), 

J B 



J 



(n+1) 



J 



n)2 



J 



(n) 



exp(-a B ), 



(44) 



(45) 



where a depends on the ratio Jx /^b - This leads to 
the modified recursion equation for as, 



where a is a function of X 



(n) 
AB 



X 



(n) 



X 



(«) 



(46) 



If the 



function a(X) is approximated by a linear function of 
X, we have 



X 



(n+l) _ 



= MX in) - (X 
«(*a°b), 



(n) 
AB 



x abH 



which can be rewritten as, 

x (»+i) = M m X {n) - a n 



(47) 



(48) 



where 



M m = 



7 = 



2 -7b -1+7s 

«(^a b)-4°b7, 
(7a,7b)- 



The iterative solution is given by, 
X (n) = M m X (n_1) 



a. 



= M n ~ 



x 1 y, u ' - 



(49) 



In this case, the equation corresponding to (17) does not 
hold for M m . By elementary algebra, we find the eigen- 
values of M m are 1 + 7b a (7ba = 7b — 7a) and unity. 
Therefore, if 7b > 7a, grows as, 

x (n) 



A"w m with A D 



1 + 7b a (50) 
This 



where it m is the eigenvector corresponding to A 
implies 



InJ 



(n) 



InJ, 



(n) 



V' 



(51) 



The effective value of 7ba, which depends on X = 
ln(JA/ Jb), is estimated in the following way. (In the 
remainder of this subsection, we omit the superscript 
("^ and distinguish the renormalizcd values by '.) Let 
us concentrate on the case Ja < Jb, because after the 
second step this inequality always holds. To determine 
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Fig. 4. The X-dependcncc of 7ba calculated from finite cluster 
numerical diagonalization. 



a a for J a < Jb, wc numerically diagonalize the cluster 
BAABABAAB. This is the largest cluster which reduces 
to a single A'-bond after decimation. Therefore the ef- 
fective coupling J' A of the A'-bond is estimated from the 
singlet-triplet energy gap of this cluster. The parameter 
a a is determined through the formula (44). To deter- 
mine «b, the same procedure is repeated with the clus- 
ter BABAABAABAB using the formula (45). By the nu- 
merical differentiation of as — a a with respect to X , the 
effective value of 7ba is obtained as shown in Fig. 4. It 
is seen from this figure that 7ba > 0. Therefore, starting 
from finite Ja/Jb, in the early stage of rcnormalization, 
the effective coupling rcnormalizcs as eq. (51) similarly 
to the case of XY model with relevant apcriodicity. 8 In 
terms of the system size dependence, this implies that 
the finite size energy gap behaves as, 



lnA£ 



InJV 
31n<£ 



A 



with 



(52) 



lnA m 



for relatively small systems. This is the reason why we 
have observed the behavior (52) in DMRG calculation in 
rcf. 9. As the renormalization precess proceeds, the ratio 
Ja/Jb decreases and 7ba tends to zero. Therefore, the 
formula (38) is recovered in the true asymptotic limit. 

4.2 XXZ case 

4.2.1 XY-like case (|A| < 1) 

For |A| < 1, the anisotropy parameters are always 
renormalized to Aa = Ab = as is obvious from the 
recursion relations (8) and (9). Therefore the XY type 
behavior AE <~ N~ z with finite dynamical exponent z is 
always expected. 

The exponent z is determined by substituting the fixed 
point values of Ja/Jb into (28) as, 

T (oo) 

(53) 



ln- B 



31n^ j(°°) 



From (6) and (7), the recursion formula for the ratio ^-f^j 




Fig. 5. The shift Az of the dynamical exponent z relative to the 
XY value zxy plotted against A. 



is given by, 



J 



(n) 



J 



(n-1) 



B 



J 



(n) 



(n-1) 



J 



-(I + Ab^- 1 ') 



(54) 



'A U A 

whose solution is given by, 



'B 



J 



(1) n-1 



J 



(n) 



J 



= TAB 



^Ila + AB 



(fch 



(55) 



fc=0 



where we have taken into account A^ = A^' = A. 
Therefore the dynamical exponent z is renormalized with 
respect to the XY value zxy as 

z = zxy + Az 
2 



,(0) 



ZXY = 



Az 



31n0 
2 

Sh4> 



lnr A j 



In jj(l + A B 



(fc)> 



(56) 



The values of Az are calculated for —1 < A < 1 as 
shown in Fig. 5. For < A < 1, the dynamical exponent 
is enhanced and diverges at A — 1, while it is reduced 
for —1 < A < 0. It should be also noted that Az does 
not depend on the ratio Ja/Jb- 

In the region, < A < 1, the Fibonacci modulation 
is always relevant within the WCRG calculation. There- 
fore, the system parameters are renormalized from the 
weak modulation regime to the strong modulation regime 
where our RSRG approach is valid. After enough steps 
of renormalization, the Fibonacci XXZ model with < 
A < 1 is renormalized to the XY chain with Fibonacci 
modulation and the energy gap scales as AE ~ N~ z with 
nonuniversal dynamical exponent z. This looks again in 
contradiction with the DMRG data in ref. 10 which are 

well fitted by \nAE . To clarify this point, we 

carried out the numerical RSRG calculation using the 
clusters used in the last section. The effective coupling 
and anisotropy parameters are determined from the ex- 
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Fig. 6. The renormalization step dependence of the energy scale 
AE obtained by the numerical renormalization group method 
with (a) Jb/Ja = 2 and (b) Jb/Ja = 3. The renormalization 
step n should be interpreted as lnAT/ (31n</>) in terms of the system 
size N. 

citation spectrum of the clusters in each step of renor- 
malization. The results are shown in Figs. 6(a,b). For 
< A < 1, the effective coupling scales as N~ z for large 
enough TV. However, for small system sizes, the behav- 
ior of AE is close to the isotropic case for non- vanishing 
A. This is clearly observed for relatively large A in Figs. 
6(a,b). This explains why the DMRG data for the energy 
gap was fitted well by the formula InAE ~ iV"\ 10 

In the region — 1 < A < 0, the situation is more 
complicated. For large enough tab, the anisotropy A is 
again renormalized to small values. Therefore, the Fi- 
bonacci XY line is approached and the scaling relation 
AE ~ N~ z with nonuniversal dynamical exponent z 
should hold. 

According to the WCRG analysis by Vidal and 
coworkers, 12 ' 13 however, the Fibonacci modulation is ir- 
relevant for —1 < A < as far as the modulation am- 
plitude is small. Therefore the ground state should be 
renormalized to the Luttinger liquid state without aperi- 
odic modulation. This looks inconsistent with our RSRG 
conclusion. However, as A approaches —1, even for large 
tab , the perturbation approach breaks down because the 
energy denominators in (3) and (4) becomes larger than 
the strength of perturbation. The limit of reliability of 
the present strong coupling approximation is given by 



Fig. 7. The TV-dependence of < ln(l/A) > of the Fibonacci XXZ 
chain with negative A for (a) Ja = 1 Jb = 2 and (b) Ja = 1 
Jb = 3 by the DMRG method plotted against IniV. 

tab > r CI ~ (1 + A)" 1 . For r AB < r cr , the WCRG 
prediction is more reliable and we expect the Luttinger 
liquid behavior AE ~ A^" 1 , C ~ T and \ ~ const. For 
tab > r CT , the marginal behavior AE ~ N~ z , C ~ T 1 / 2 
and x ~ T 1/z_1 is expected. Actually the DMRG data 
for Jb/Ja = 2 and 3 with negative A shown in Fig. 7 
confirms the clear deviation from the Luttinger liquid 
behavior AE ~ A^ -1 for small negative values of A. The 
details of the DMRG calculation are as described in refs. 
9 and 10 except that the maximum number of states re- 
tained in each subblock is 160. This value is required due 
to the large modulation strength. 

4.2.2 Ising-like case : (|A| > 1) 

For A < —1, it is obvious that the ground state is 
the ferromagnetic state. For A > 1, the anisotropy Aa 
and Ab are renormalized to infinity as seen from (8) and 
(9). Therefore the ground state is the Neel state. How- 
ever, in this case, the perturbational recursion equation 
breaks down because J^ A^ becomes larger than J^ 
for large enough n if our recursion equation is naively 
solved. This is reasonable because the Neel state is not 
adiabatically connnected to the dimer state which is the 
starting point of our RSRG scheme. Therefore, it is not 
possible to make further quantitative predictions in this 
regime within the RSRG scheme. 
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Fig. 8. The schematic ground state phase diagram of the Fi- 
bonacci XXZ chain. 



4-2.3 Phase Diagram 

All the above arguments are summarized as the 
schematic ground state phase diagram in Fig. 8. The thin 
solid curve corresponds to the line tab = r„ ~ (1 + A) -1 
which separates the Luttinger liquid phase and the Fi- 
bonacci XY phase. The thick solid line is the isotropic 
Fibonacci Heisenberg phase with scaling AE ~ N~ KlnN . 
The thick broken line belongs to the Luttinger liquid 
phase. The thick dotted line is the fixed line of the 
marginal Fibonacci XY state. The Neel state and the 
ferromagnetic state appear for A > 1 and A < — 1, re- 
spectively. 

5. Summary and Discussion 

The ground state properties of the S = 1/2 Fi- 
bonacci XXZ chain arc investigated by means of the 
RSRG method. In the isotropic case A = 1, it is found 
that the energy gap scales as ]\[- KlnN which implies the 
logarithmic divergence of the dynamical exponent. For 
— 1 < A < 1, the RSRG flow converges to the Fibonacci 
XY-line. Therefore, the marginal behavior with finite 
dynamical exponent z is expected as far as the RSRG 
scheme is reliable. With the help of the WCRG results 
by Vidal and coworkers, 12, 13 which is valid in the weak 
modulation regime with — 1 < A < 1, the ground state 
phase diagram is determined. The previous DMRG re- 
sults by the present author arc reinterpreted from the 
viewpoint of the present theory and understood as the 
finite size crossover effect. 

We have found a new quantum dynamical critical be- 
havior (38) which was so far unknown in the field of 
quantum many body problem. Similar 'singular dynamic 
scaling' is, however, known since 80's for the classical 
Ising model on the percolation clusters with Glauber 
dynamics. 20-22 In spite of the geometrical self-similarity 
common to this classical model and our quantum model, 
they look very different in many aspects. Although the 
underlying physics which leads to the similarity in the 
dynamical critical behavior of both systems is unclear, it 
would be possible that a kind of hidden symmetry under- 
lies in the problem of the Ising model on the percolation 



ABAABABAABAABABAABABAABAABABAABAAB * * * * 

4 4 4 4 4 4 44 

Fig. A-l. Construction of the Fibonacci sequence from A3 and 



clusters considering that the singular scaling (38) in our 
model only occurs in the isotropic case (SU(2) symmetric 
case). Further investigation on this point might lead to 
a more profound understanding of the dynamical critical 
behavior of both systems. 

It is well known that the RSRG method is suitable for 
the random spin chain problem. 19 ' 23-25 The present cal- 
culation demonstrates that the RSRG method is also a 
powerful tool for the investigation of the quasiperiodic 
chains. Actually, after this work is almost completed, 
the preprint by Vieira appeared in e-print archive in 
which some of the present results are derived. In addition, 
Vieira satisfactorily applied this method to the Heisen- 
berg chains with relevant aperiodicity. Similar approach 
is also applied to the two dimensional quasicrystal. 18 We 
thus expect the RSRG method is widely applicable to 
various problems in the field of quantum magnetism in 
quasiperiodic systems. 

The author would like to thank Dr. C. L. Henley for 
drawing his attention to refs. 20-22 and for enlightening 
comments. The numerical computation in this work has 
been carried out using the facilities of the Supercomputer 
Center, Institute for Solid State Physics, University of 
Tokyo and the Information Processing Center, Saitama 
University. This work is supported by a Grant-in-Aid 
for Scientific Research from the Ministry of Education, 
Culture, Sports, Science and Technology, Japan. 

Appendix 

In this appendix, we prove that our decimation pro- 
cedure corresponds to the deflation process of the Fi- 
bonacci sequence. 

Let us build a Fibonacci sequence from the third gen- 
eration building blocks A 3 = ABAAB and B 3 = ABA 
as shown in Fig. A-l. The thin upward arrows indicate 
the boundary of these building blocks. Following the dec- 
imation process described in the text, the spins survive 
in the middle of the sequence AA as indicated by thick 
upward arrows. 

From the construction rule, the ABA sequences which 
are common to the left ends of the A 3 and B 3 arise at 
the right side of every block boundary. If the next block 
is A3, the third position to the right of the boundary is 
always AA. If the next block is B 3 , the next to next block 
is always A 3 because there appear no successive B's in 
the Fibonacci sequence. Therefore in this case also the 
third position to the right of the boundary is AA (next 
boundary). This implies that there always appear the 
sequence AA at the third position to the right of the 
boundary. Because these are the only cases in which AA 
sequence is allowed, at the third position to the left of 
the AA-sequence, we have always the block boundary. 
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Thus, the sequence built by the substitution rule A 
— > ABABA and B — ► ABA, which is the inverse of our 
decimation procedure, is a Fibonacci sequence lacking 
in the first three letters ABA. Therefore, if we apply our 
decimation procedure to the Fibonacci sequence, we have 
again a Fibonacci sequence with an additional B' at the 
leftmost position. 
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